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ABSTRACT 

The NASA/mSFC - Multilevel Diffusion Model (MDM) Version 5 used ty 
the Environmental Effects Office at JSC was modified to Include features of 
more recent versions. The MDM was used to predict In-cloud HCl concentra- 
tions for the April 12 launch of the Space Shuttle (ST.^ l). The maximum 
centerline predictions were compared with measurements of maxlmvun gaseous 
HCl obtained from aircraft passes through two segments of the fragmented 
Shuttle ground cloud. The model over-predicted the maximum values for 
gaseous HCl in the lower clou.'* »egment and portrayed the sane rate of 
decay with time as the observed values. However, the decay with time of 
HCl maxlmums predicted by the MDM was more rapid than the observed decay 
for the higher cloud segment, causing the model to under- predict concentra- 
tions which were measured late in the life of the cloud. The causes of 
the tendency for the MDM to be conservative in over-estimating the 
HCl concentrations In the one case while tending to under-predict concentra- 
tions in the other case are discussed. 

Further comparisons of the MDM predictions for in-cloud HCl concen- 
trations were made for Titan III launches in which aircraft measurements of 
HCl were available. These comparisons indicated that the model is conserva- 
tive and over-predlcts the maximum HCl concentrations early in the cloud's 
his''.ory. Results for in-cloud HCl concentrations for some of the meteorologies 
characteristic of Cape Canaveral axe presented. 
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IMTRODUGTION 

The prlmairy objectives of the work reported on here were toi 
(l) Develop the capabilities of the NASA/MSFC Multilayer Diffusion Model 
(MDM) Version 5 to obtain In-cloud predictions of HCl concentrations for 
the Space Shuttle, and (2) To use the MDM for comparisons with the field 
measurements on the first Shuttle launch for the purpose of model valida- 
tion. 

Additional objectives of the woric reported on here Included 
obtaining Information ont (l) The capabilities of the MDM for predicting 
In-cloud concentrations of HCl for Titan III launches, and (2) The In- 
cloud HCl concentrations which may be expected to be encountered for Space 
Shuttle launches at Cape Canaveiral under average meteorological conditions. 

The results reported here include the procedures used to implement 
the In-cloud prediction capabilities for the MDM on a PDP-ll/45 at the 
Environmental Effects Office at JSC. The documentation Is given for the 
changes made In the MDM which allow the selection of Version 6 Shuttle 
parameters already in the program. Version 7 parameters were added to 
the selection of options available for Space Shuttle launch parameters. 

A comparison of the HCl predictions by the MDM for a given meteorology 
using the three versions of Shuttle launch parameters indicated that there 
was an Insignificant difference between their predictions. 

The MDM validation for the April 12, 1981, Space Shuttle launch 
(STS-l) showed that the model predictions for the lower cloud segment 
in an unstable setting closely portrayed the decay of HCl with time but 
over- predicted the magnitude of gaseous HCl. One factor which could be 
expected to cause the measured HCl concentration to be below model 
predictions is that this segment of the cloud had a high relative humidity. 
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In-cloud neasurements indicated that this resulted in a large fraction of 
HCl being involved in aerosol rather than gaseous form. The MDM predictions 
for the higher cloud segment which drifted westward from the launch site 
were less than the meaisured values. The upper cloud was not sampled until 
eifter passes through the lower cloud had been completed or about 50 
minutes after launch. The MDM predictions for HCl which were for a much 
later period in the cloud's history were lower them those observed. This 
may be due in part to the assumption in the MDM that the diffusion processes 
continue at a constant rate throughout the cloud history. The diffusion 
rate is determined by the standard deviation in the horizontal wind from 
values near the surface. As the laiuich cloud enters a more stable 
environment, as was the case for the upper cloud segment, this assumption 
would tend to cause the MDM to over-estimate the rate of decay of HCl within 
the cloud. Measurements also indicated that essentially all the HCl was 
in gaseous form in this cloud which had low relative humidity. 

Predictions for In-cloud HCl concentrations for Titan launches 
indicated a tendency for the MDM to over-estimate the concentrations. The 
in- cloud HCl concentrations for Shuttle launches predicted by the MDM for 
the standard meteorologies at Cape Canaveral closely parallel those for 
Titaui launches for the same atmospheric conditions. 
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MODIFICATIONS OF THE MDM 

The Environmental Effects Office at JSC has used the NASA/MSFC 
Multilayer Diffusion Model Version 5 137 Dumbauld and Bjorklund (l) to 
predict surface concentrations of Hd for rockets using solid fuel boosters 
(Glasser et al., 2, 3)> In Its original form the MDM was run on UNIVAC but 
it was modified In 1977 by Joe Yoder to run on a PDF 11/45. Further 
changes were made In 1979 by larry Ray of JSC Including the addition of 
Version 6 parameters (Dumbauld and Bjozklund, 4) which Included new 
values for the Space Shuttle. 

Attempts to use the MIM In its earlier and modified forms at JSC 
to give HCl concentrations at levels other than the surface were not 
successful. To Implement this capability of the MDM, it was necessary to 
modify two statements in subroutine SETUP DAT., l.e., NPTS which 
Identifies the number of levels in the cloud at which concentrations 
are desired and ZZL which is the parameter for the height of these levels 
in meters. Invisible errors In thiB routine, l.e., blank spaces which did 
not show In printouts frustrated earlier attempts to obtain in-cloud con- 
centrations. With these modifications the MDM will output HCl concentra- 
tions at any level up to the cloud stabilization height. 

Because of the great length and complexity of the MDM, even small 
modifications can be difficult to accomplish or can cause unexpected 
problems. This is particularly true if more than one segment of the 
program is involved. Richard Roenfeldt made changes in the MDM which 
implemented the Version 6 constants. For example, even though Version 6 
constants had been put into the program and when called, the program 
indicated they haul been used, careful checking on the outputs, however, 
showed Version 5 constants were used in all cases. Appendix A gives 
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the documentation on changes required to implement Version 6 and Version 7 
constants. Richard Roenfeldt msde these and other changes in the MDM 
relative to this report. Table 1 gives the most recent constants for 
the MDM which have been used in the REEDM version of the MDM being used 
at KSC. These constants do not differ substantially from Versions 5 and 
6 constants £uid did not cause significant differences in HCl predictions 
when run on identical cases. 
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MDM PREDICTIONS FOR 3TS-1 

The most important aspect of the work reported here is the valida- 
tion of the MDM for predicting in-cloud KCl cozicentrations for Space 
Shuttle launches. The results of field measur«ients of HCl from air- 
craft flights through the ground launch cloud for STS-1 as reported hy 
Sebacher et al. (5) were used in this validation. 

The Data 

The meteorological sounding at launch time (T-O) for STS-1 on 
April 12, 1981, is given in Table 2. This data was obtained by telephone 
transmission through the interface with the KSC computer and was used in 
the MDM predictions reported here. An abbreviated version of the 
meteorology giving a few of the levels is provided in Table 3- A graph 
of the temperature and dewpoint tonperatures as a function of height from 
these tables is shown in Figure 1 . On this same Figure the height predicted 
bythe MI^ for stabilization of the laiinch cloud (1187 m) is given. The 
flights through the fragmented launch clouds A which ranged from 850 m to 
900 m, and for Cloud B from 1600 m to 18?0 m are also given for reference 
purposes . 

The temperature sounding of Figure 1 clearly shows a shallow 
surface inversion and a moderate upper level Inversion and stable layer 
extending from 3256 feet to 7000 feet. This type of sounding Is ch 2 iracterlstic 
of weather regimens for the Cape in which the Bermuda High extends over 
the Florida Penisula. Subsidence In the high pressure area produces the 
Inversion and stable layer at upper levels. This stable layer Is 
responsible for suppressing the observed stablizatlon height for Cloud B 
but the Inversion is not Intense enough to suppress the launch cloud to 
the level predicted by the MDM. 
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Ebcaaination of tho wind directions In Table 3 Indloate a vertical 
shear of the horizontal wind throughout the mixing layer and across the 
inversion. This shear helped to fragment the Space Shuttle launch cloud 
which was observed to stabilize in five segments each at a different 
height . 

Figure 2 provides a '/lew of this wind shear at the launch site. 
Arrows represent the magnitude and direction of wind for the heights In 
Table 3« The wind direction near the 1000 foot level is generally north- 
ward nearly along the shoreline of the Cape and corresponds to the direction 
the lower cloud ( Cloud A) was observed to travel. Figure 3 gives a rough 
sketch of the path taiken by Cloud A for the first 24 minutes of observa- 
tion. The wind direction at the 6000 foot level is towaird the west and 
more nearly corresponds to the direction that the upper cloud (Cloud B ) 
was observed to travel. The direction of cloud movement predicted by the 
NDH is intermediate to the two clouds and is represented ty crosshatching 
in Figure 2. 

The fragments of the Shuttle launch cloud were observed to reach 
stabilization height 8 minutes after launch. Sampling of Cloud A for 
HCl gases and aerosols and for parti collates begin at 8.6 minutes after 
launch for Cloud A and continued at 2 to 5 minute intervals until 45 
minutes after launch. The higher cloud waus similarly sampled from 49 
minutes until 2 hours and 8 minutes after launch. Ebcamples of the HCl 
mesisurements for aircraft paisses through the upper and lower clouds are 
given in Figure 4 (from Sebacher et al., 5)* ^'be low altitude segment, 

Part (a) of Figure 4 shows that Cloud A has a high relative humidity 
and has a small fraction of HCl in gaseous form while much of the HCl is 
contained in aerosol forji. The high altitude segment, Part (b) of Figure 4 
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•hows that Cloud B has a low ralatlve humidity and has nearly all of 
Its HCl In non-aerosol form. The aaxlauffl value of gaseous HCl for each 
pass through Clouds A and B Is given In Table 4. The maximum total HCl 
are plotted In Figure 5 (Sebacher, 5) along with particulate concentra- 
tion, relative humidity aiul twperature. It should be noted that these 
values are plotted as a function of time after launch. This adds an 
element of uncertainty when making comparisons with HCl predictions 
from the NDN since thMe predictions are given ais a function of distance 
from the point of launch. 

The in-cloud HCl predictions computed here were obtained from 
the HDM using the meteorology from Table 2. In-cloud concentrations 
were computed for four different levels corresponding roughly to the 
upper and lower limits of adrcraft sampling heights for Clouds A and B 
(see Figure 1). Values of maximum centerline HCl for 850 m and 900 m 
were obtained for the lower cloud and for 1600 m and 1800 m for the upper 
cloud (Table 5)* 

The maximum peak HCl predictions for the lower cloud at the 85 O 
and 900 meter levels given In Table 5 differ by less than 1^, while those 
for the upper cloud differ by less than 10^. The maximum peak (centerline) 
HCl concentrations from Table 5 ue plotted In Figure 5 . HCl measurements 
recorded in Table 4 were madt' as a function of time in reference to the 
launch. The MIBl predictions, however, axe output as a function of 
distance of the launch cloud _'rom the launch site. In order to make a 
comparison of these HIM predictions with the HCl measurements, it is 
necessary to make some assumption rel>*tive to the equivalence between 
the time from launch and distance of the launch clo\id from the launch 
site. The most reasonable aissumption would be to consider that the 
cloud fragments move with a speed equal to the average wind speed of 
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th« layer at which thw particular cloud stAbl.llsad. From Table Z the 
wind speed decreactes from 12 knots at 2000 feet to 9 knots at the 3000 foot 
level. Since Cloud A drifted northward at allttudes fron 650 m«>ters 
(2(133 feet) up to 950 meters (3fH7 feet). It would be reasonable to 
assume that It e?;perlenced an average on the order of 10.5 knots 
(5>9’ m/sec). The second cloud segment was observed to drift westward at 
altitudes from 1350 moters (4,429 feet) up to 1380 .leters (6,168 feet). 

From Table 2 the wind increases from 8 knots at the 4000 foot level to 
16 knots at the 6551 foot level making it reasonable to assume an average 
wind on the order of 12 knots (6.17 m/sec) for Cloud B. It is necessary 
to add to the values output by the MOM the amoiint of time elapsed from 
lauinch to cloud stabilization which was at 1250 m and 2500 n downwind from 
the launch site according to the MOM. Using an average wind speed for the 
rising larnich cloud of 10.5 knots given a time to cloud stabilization for 
the lower cloud of 5 mimites, 21 seconds and a time of 7 minutes, 43 
seconds for the upper cloud. This is close to the 8 minutes to cloud 
stabilization that was reported to be observed by Sebacher (5)« It will 
be noted later in this report that a shift of the time scale by several 
minutes In either direction will not significantly alter the conclusions 
reached relative to the comparison of observed and predicted HCl concentra- 
tions. The 'ralues for the correspondence between time and distance 
scales for Clouds A and B using the assiinptions discussed are tabulated 
in Table 6 along with HCl predictions from Tabic 5 aukl are used in 
Figures 7 and 8 to compare HCl observations and predictioiis. 

The Analysis 

In Figure 7 the MIM predictions for peak centerline HCl concentra- 
tions given by the solid line exceed the peak values of gaseous HCl 
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represent ed by squares which were through the lower cloud. The agreement 
between the meignltude of the observed and measured values of gaseous 
HCl is fair considering the uncertainties inherent in both methods of 
determinirj it. 

The rate of decay of HCl with time is in p;^xticularly good agree- 
ment for both predicted and measured values. The lower cloud is in a 
region where the atmosphere is less stable than where the upper cloud is 
located. This may be determined by looking at the temperature profile in 
the plot of the MET sounding in Figure 1. The rate of decay of HCl 
concentration as determined by the MDM is largely a function of the 
standard deviation of the horizontal wind speed ( ) as used in the 

diffusion c^J.culations . This parameter was obtained from the KSC computer 
which calculates <r using atn objective routine that analyzes the variances 
in wind direction. The value of <T = 13 which was used is relatively 
large as parametric studies (Glasser, l) have shown. This value of ^ 
would appear to be representative for HCl concentration decay in the region 
below the upper level inversion shown in Figure 7. 

The magnitude of HCl concentrations predicted by the MDM has been 
shown to be conservative in other studies which have used it to predict 
surface concentrations of HCl for Titan launches. The over-prediction of 
in-cloud HCl would also be expected because of conservative assiunptions 
which have been built into the MDM. Another factor which would tend to 
cause the predicted HCl values to be larger thaui the measured values 
for this particular case is the large amount of HCl that is in the 
aerosol form. The measured values for total HCl (gaseous plus aerosol) 
is given in Table 4 and plotted in Figures 4, 5, and 7. In Figure 7 it 
can be seen that the MDM predicted value lies about midway between 
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gaaeous and total HCl conoentrationa. The rate of decay of total HGl 

closely parallels the rate of decay of the predicted and measured HCl in 

gaseous form. 

The agreement between measured and observed HCl could also be 
affected by the assumptions used In relating measured sampling time to 
predict cloud position. The data on observed cloud position given In 
Figure 3 indicates Cloud A Is 5 ^ troa the launch site In 10 minutes 
and 10 km from the pad 39 In 24 minutes. This compares with the calcula- 
tions used on the MDM predictions given In Table 6 whore at 5 i™ the time 

Is 13 minutes, 36 seconds and at 10 km the time Is 30 minutes, 32 seconds. 

To bring the time assumed for the MDM cloud position Into agreement 
with the cloud positions in Figure 3 would require a subtraction of 
about 6 minutes which would have the effect of shifting the MDM prediction 
to the position of the dashed line In Figure ?• This is not enough to 
affect the analysis of the comparison of observed ajid measured HCl 
values given here. 

In Figure 8 the MDM predictions for HCl In Cloud B, represented 
by a solid line, are compared to measurements of gaseous and total HCl 
concentrations. The MM predictions in contrast to those for the lower 
cloud significantly under- predict by a factor of about 3 the gaseous HCl. 
The measurements of gaseous and total HCl also do not display the decay 
with time predicted by the model. In fact, the gaseous HGl values decay 
relatively slowly over the 70 minutes of sampling time as Indicated by 
the dashed line In Figure 8. 

The reasons for the lack of agreement are probably related to the 
fact that the upper cloud has entered a stable environment above the 
Inversion (note Figure 1). In this environment mixing processes are 
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Inhlblted while the lOf essentially assumes the same rate of decay 
established by the choice of ft* In the surface environment. It would have 
been useful to have HCl concentration measurements of Cloud B early In 
Its history to check on the role of the decay rate In this over-prediction 
by the NDN. 

Another difference between the lower <ind upper clouds is that Cloud B 
had a low relative humidity causing the HCl concentration to be almost 
entirely In the gaseous phaae. The measurements of total and gaseous 
HCl plotted In Figure 7 show a great degree of variability perhaps 
suggestive of the difficulty In making accurate measures under these 
circumstances. The error range In these measurements was provided by 
Richard Bendura of LRC aa i 20^ with a precision of measurement of 0.5 PPm. 
The variability of the data could adso be related to the difficulty of 
aircraft sampling when the cloud has become diffuse with the passage of 
so much time. 

One problem with the use of the MDM for making these predictions 
is certain to cause the HCl values to be under-estimated is the following. 

The MDM will not compute HCl concentrations above the mixing height which 
must be chosen subjectively prior to running the program. From Figure 1 
the height of the surface mixing layer is clearly at the base oi the 
upper level Inversion. In order to have the MDM calculate concentrations 
above this level, it was necessary to assume the mixing would occur 
throughout the layer from Cloud B to the surface. This assumption is 
not realistic and causes the concentrations of HCl to be reduced at every 
level. It is, therefore, quite probable that the under- predict ion of 
HCl concentrations in the upper cloud axe related to problems inherent in 
the MDM which prohibit It from more realistic modeling changes encountered 
in the real atmosphere. 
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ADDITIONAL COMPUTATIONS 

The proceeding material completes the report as required for the 
original grant application (Appendix B) . During the period of work on 
Phase I| which was for one month at JSC, the objectives of Phase II, 
which was projected to require six months, were also primarily completed. 
Because the capabilities of the MDM for predicting In-cloud HCL concen- 
trations were Implemented so quickly, additional In-cloud data was 
developed while awaiting the Space Shuttle laxmch and while work on 
modifying the MDM proceeded. Since that data does not appear to contribute 
much that Is new In the way of Insights Into the subject of HCl concen- 
tration predictions amd since complete documentation of these data would 
produce a very lengthly report of potentially little merit or interest, 
the nature of the data will only be briefly summarized to Indicate what 
is available. 

One area of Interest was the in-cloud HCl concentration predictions 
for rocket launches using solid fuel boosters for which in-cloud HCL measure- 
ments had been made. 

In-cloud predictions were made for Titan III launches for a 
number of different weather regimens including one using MET data for the 
May 20, I975i Tltam launch for which HCl measurements were available. 
Examination of the case appeared only to confirm the results by Rudolph 
(6) which indicated the tendency of the MDM to over-predlct HCl values 
for Titan and Delta launches for the period studies, 1973-1978. An 
example or the type of data developed for Titan launches is given in 
Table 7 for February 27, 1965, which gives the peak HCl concentrations 
for each 100 meter level through the mixing layer. The highest HCl 
values center on 700 meters while the cloud stabilized at 900 meters. 
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The highest HCl values also occur closest to the launch site, 7500 meters, 
while at levels away from cloud center the peak HCl concentration occur 
at greater distances. Indicating the time lag as diffusion takes place both 
upwaixi and downward. 

Another area of Interest Is the effect of the different weather 
regimens found to occur at Gape Canaveral as presented by Slier (?), on 
In-cloud HCl concentrations for Shuttle launches. In this analysis, the 
weather regimens In which the subtropical ridge dominates the weather of 
the Cape are classed as Al, A2 or A3 depending on whether the pressure 
center lies close, south or north of Canaveral. The April 12 Shuttle 
launch occurred under Al conditions which have the greatest probability 
of occurring (20^) of any of the six weather types characteristic of the 
region. Under these conditions the probability of onshore tiansport of 
the launch cloud Is over 90^. A vertical profile of Al weather from 
the day of the Shuttle launch, April 12, 1981, Is given In Table 8. 

The vertical profiles of peak HCl from MDM predictions In Table 8 
caul be related to the results of the Shuttle launch In the previous 
discussions. In the MDM predictions with the 3000 foot mixing layer 
assumption, the HCl Increases to the 1300 m level to 65 Ppm while In the 
caise with the 3750 foot mixing level, the HCl Increases to 26 ppm at the 
600 m level. The effect of changing ^ from 4.5 to 9*0 Is to markedly 
decrease HCl concentrations at the cloud center and to Increase values 
below the doud Indicating rapid mixing of HCl throughout the layer. 

This demonstrates the effect of an Increase In ^ on decreasing the HCl 
predictions by the MDM. The over-prediction of HCl for the upper cloud 
In the Shuttle launch has been considered to possibly be attributed to 
the or of 12.0 used In those predictions. 
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MDH predictions for HCl concentrations were made for several 
levels In the vertical for MET conditions representing each of the 
weather types developed by Slier (?) . These predictions were made for 
both Shuttle aind Titan launch parameters. Although this rather large 
amount of data represents a kind of climatology of In-cloud HCl concen- 
tration predictions. Its value Is somewhat reduced by the lack In 
uniformity In sissumptlons . This Is because the data was developed over 
several months while changes were being made In the MDM Itself. 
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OONCmSIONS 

ThiB work represents a first attempt to compare In- cloud HCl con- 
centration predictions to In-cloud aircraft measurements of HCl for the 
Space Shuttle launch. The inadequacy of the NASA/MSFC MDM to accurately 
portray the actual complexities of the diffusion process and particularly 
to cope with the effect of changing conditions which rocket launch clouds 
encounter as they drift from the site are well known and have been given 
consideration In numerous studies. If there Is a general conclusion from 
the work presented here, It Is that In spite of the numerous experimental 
and theoretical difficulties In obtaining the In-cloud HCl concentrations, 
the agreement Is at least within an order of magnitude. 

The fragmentation of the Shuttle launch cloud on the April 12, 

1981, launch presents a serious difficulty for the Mm at the onset since 
only simple cloud geometries are assumed. In spite of these difficulties, 
the decay rate of peak HCl concentrations In the lower cloud ao:e well 
portrayed by the MM! and are only slightly over-predicted. The over- 
predictions may be understandable as discussed because of the significant 
amount of HCl which Is in aerosol form due to the high relative humidity 
of the lower cloud. 

The decay of HCl concentrations predicted by the MDM for the 
upper cloud is much more rapid than observed over the ?0 minute sampling 
period. As discussed, this could be related to the use of a standard 
deviation of the horizontal wind direction ( <T") that is appropriate for 
estimating the diffusion processes In the lower cloud which is in an 
unstable environment. The upper cloud, however, is in a region of 
generally high stability which reduces mixing. This could adso account 
for the magnitudes of Hd being under-predicted particularly since the 
upper cloud was not sampled until about 50 minutes had elapsed. In 
general, it is apparent from this study that the MM! can produce in- 
cloud HCl values that fall within a reasonable range of measurement. 
Comparisons of MDM HCl concentrations with surface HCl measurements show 
less agreement since studies Indicate it over-predlcts by an order of 
magnitude or more. 

If more refinement is required in the knowledge of in-cloud 
HCl, it is likely that both the model and the measurements will have 
to be improved. 
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TABI£ 3* Abbravlatad Radlosonda Data for STS-1 

Data for Cape Caiiaveral on 1212 Z April 12, 1981, coxreepondlng 
to the launch of STS-1. Data was provided by Richard Bendura of litC but 
corresponds to excerpts f^rom the ccaplete MET data set In Table 2. 


A1 titude, 
ft. 

Direction, 

Wind 

^ ^Speei^ 

Temperature 

*C 

Dew Prtt. 

«c 

Pressure 

mb 

Rel . 
Hum. , % 

16 

110 

4 ' 

, 17.0 

15.9 

1023.4 

93 

1000 

136 

12 

19.1 

14.8 

988.46 

76 

2000 

142 

12- 

16.2 

13.5 

953.98 

34 

3000 

136 

9 

14.1 

11.7 

920.41 

86 

4000 

099 

8' 

15.1 

- .6 

887.90 

37 

5000 

079 

12 

15.3 

-2.2 

856.56 

30 

6000 

074 

15 ' 

14.2 

1 

ro 

• 

826.26 

31 
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TABU 4. In-doud HCl Data for STA-1 

The HCl data obtained by aircraft paesee through the Shuttle 
Ivcroh cloud fragaenta A and B. The values reinresent the peak HCL eoncentra 
tiona from each pass aa provided by Richard Bendura of URC, and used in 
reports by Sebacher et al. (3). 

Tiae froa Total HCl^ 

Cloud Aircraift Launch Aerosol Gaseous HCl 

Fkagaent Pass No. (aim sec) Gaseous (ppa) (ppn) 


1 

9t00 

17.5 

3.6 


2 

14i39 

11.5 

2.2 

A 

3 

19i32 

5.5 

1.4 

Lower 

4* 

23i33 

0.8 

0.2 

Cloud 

5 

27i33 

4.2 

1.0 


6 

31i57 

5.0 

1.6 


7 

'>i48 

3.0 

0.6 


8 

-9i18 

3.5 

1.2 


9 

43i25 

3.2 

1.0 


10 

54i26 

6.5 

4.6 


11 

57i50 

6.5 

3.8 


12 

62i40 

3.5 

3.8 


13 

69i47 

5.2 

4.5 

B 

14 

84i59 

2.5 

2.5 

Upper 

15 

91i58 

3.0 

2.1 

Cloud 

16 

93«41 

2.6 

2.6 


17 

95«16 

3.0 

2.4 


18 

99«55 

1.6 

2.0 


19 

I05il4 

1.8 

2.0 


20 

10915 '* 

4.2 

3.6 


21 

114i28 

2.3 

2.6 


22 

II 81 50 

2.0 

2.5 


23 

123135 

1.5 

2.2 


24 

128t24 

2.7 

3.0 

* Pass 

4 was below the 

visible cloud. 




♦* KCl '/alues aire ^ 20 % or 0,5 ?I® - whichever is greater. 
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TABU! 5* In-CLoud Hd Prtdictiona for STS-1 

Th« peak uxiBum (oentorlina) Hd oonoentzmtlons fro* MDM pre- 
dictions using the MSI dsts of Table 2. These in-cloud eonoentrations 
are for the 850 and 900 seter levels for doud A and the 1600 and 1800 
aeter leiisla for doud B. 


Range - 
Distance 

Maxlaua Peak (Cmterline) Hd 
Concentration (non) at Levels 


Fsraa lAunch 
(meters) 

Q 2 SL & 

900 a 

1600 a 

1800 a 

1250 

24.00 

23.97 



2500 

14.15 

14.37 

56.53 

61. OC 

375c 

8.03 

8.11 

27,48 

29.58 

5000 

5.01 

5.03 

15.69 

16.40 

6250 

3.51 

3.52 

9.74 

10.05 

7500 

2.71 

2.71 

6.40 

6.55 

8750 

2.21 

2.21 

4.41 

4.49 

10000 

1.85 

1.85 

3.19 

3.23 

11250 

1.57 

1.57 

2.41 

2.43 

12500 

1.34 

1.34 

1.89 

1.90 

13750 

1.16 

1.16 

1.53 

1.54 

15000 

1.01 

1.01 

1.28 

1.26 

I625O 

0.89 

0.89 

1.08 

1.08 

17500 

0.79 

0.79 

0.93 

0.93 

19750 

0.70 

0.70 

0.81 

0.81 

20000 

0.63 

0.63 

0.71 

0.71 

25000 



0.46 

0.46 

30000 



0.32 

0.32 

35000 



0.23 

0.23 

1-0000 



0.18 

0.18 

45000 



0.14 

0.14 


TABUS 6. Tl»e-Dlstance and HCl Predictions for the STS-1 launch CGuOuda A and B 
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TABI£ 7. Sxaffiple Vertical Profiles of HCl from the MI^ for Titan III 
February 27 > 1965t launch 

The HGl concentrations are peak values predicted by the MDM for 
each 100 meter level. The MET data is for 051 7Z February 27 • 1965i and 
the depth of the mixing layer is 3529 feet, the cloud rise is 904 meters 
and the standard deviation in the hirozintal wind direction is 7> 

Altitude Distance From MDM Prediction 

(m) Launch (m) of Peak Hd (ppm) 

0 

11250 

0.95 

100 

10000 

1.08 

200 

8750 

1.61 

300 

8750 

3.03 

400 

7500 

8.20 

500 

7500 

16.78 

600 

7500 

25.81 

700 

7500 

29.80 

800 

8750 

9.40 

900 

8750 

8.30 

1000 

8750 

7.26 







TABLE 8. Example Vertical Profiles of HCl Ttom the for Shuttle Launch 
April I2t 1981, Launch 


The Hd concentrations are peak values predicted by the HDM for 
each 100 meter level. The MET data Is for OOOZ April 12, 1981, the day of 
the Shuttle launch. The depth of the nixing layer Is 5000 fsc't and 3750 
feet, the cloud rise 120? meters and 78? meters and the standard deviation 
of the horizontal Hind direction was 4,5 and 9> 


Al.tltude 

(m) 

Peak HCl 
MDN (ppm) 

r » 4.5 

Distance From 
Launch (m) 

Peak HCl 
MDM (ppn) 

r - 9.0 

Distance From 
Launch (m) 

0 

12500 

1.24 

5000 

3.56 

100 

12500 

1.28 

5000 

3.79 

200 

10000 

1.40 

2500 

4.48 

300 

1250 

2.88 

1250 

6.39 

400 

500 

5.34 

500 

10.57 

500 

500 

5.3^ 

500 

18.89 

600 

500 

10.67 

500 

26.44 

700 

500 

21.24 

1250 

19.77 

800 

500 

21.24 

1250 

19.01 

900 

500 

31.57 

1250 

17.13 

1000 

1250 

38.15 

Mixing height was reduced 

1100 

1250 

38.64 

from 5000 to 3750 

1200 

1250 

49.85 

feet and 

cloud rise reduced 

1300 

2500 

65.00 

from 1207 m to 787 m. 

1400 

2500 

64.00 
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PIGURE 3. Cloud Track for Qoud A of STS-1 

A rough sketch of the movement of the lower cloud fragment A J I 

from the Shuttle launch at 1212 Z April 12| 1981. The movementt which ! j 

roughly parallels the coastline of Cape Canaveral, is indicated «is a function 
of time. Sketch provided by Richard Bendura of IRC. 
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FIGURE 4. SxaaplM of In-Cloud Maaaurenents for STS-1 

Typical moasviroments for total HCl, gaseous HCl, particulate 
concentration* relative humidity, and temperature for Aircraft Pass 2 
through Cloud A and Aircraft Pass 11 through Cloud B. Graphs are from 
Sebacher et al. (3)> 
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FIGURS 5* ValuM for In-Qoud SaapUng of STS-1 

The values plotted here represent the peak measurements of total 
HCI, gaseous HCI, particulate concentration, relative humidity, and 
temperature for each pass through the upper and lower Space Shuttle 
Ground Cloud vs. the time after launch. The gratis are from Sebacher et al. 
( 5 ). 
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FIGUBS 6i MDM Pradictlons of Hd for STS-1 as a Function of Distance 

The MDM predictions of the upper fragment of the Space Shuttle 
launch cloud for the 1600 meter level Is given by the upper curve. The 
lower curve Is the prediction for the 900 meter level of the lower cloud 
fragment using the MBIT data from Table 2. 
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FIQURE 7* Measurad and Fradicted In-doud HCl Goneantratlona for 
Cloud A I STS-1 


Tha solid Una raprasanta tha in-cloud Hd ooncantratlona pra- 
dictad tha MM for tha 850 aatar laval. Tha data points oarkad with 
an X ara for total Hd Including gaaaous and aarosol. Tha data points 
oarkad with a sqxuura ara for tha aaasuranents of gasaous HCl only. Tha 
nuobars by tha data points Indloata tha flight pass numbar. *nia data 
valuas an for tha lowar cloud (a) tskan tram Tabla 4, Sebachar at al. 
( 3 ) • Tha dashad Una nprasants an adjustoant of Mm pradlctlons taking 
into account obsarvad movem«its of doud A glvan In Flgura 3« 
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original PAGE IS 
OF POOR QUALITY 

FIGURS 8. Haastirod anl Predict ad In-Cloud HCl Concentrations for Cloud B, 
STS-1 

The solid line represents the in-cloud HGl concentrations predicted 
hy the MDH for the 1600 or 1800 meter level. The data points masked xith 
an X are for total Hd including gaseous and aerosol. The data points 
nariced >dth a square are for the meMureaents of gaseous HCl only. Some 
of the data points have corresponding fli^t pass numbers adjacent to them. 
The data values are for the upper Cloud (B) taken from Table 4, Sebacher 
et ed. (5). 
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APPENDIX A 


HASA DocumenCatlon 
A. Taskrbuildlzig. thA.QJCOD modulft: 

Th« major problem with making dlanges^ to Che programs were to Cask*build 
Che system, after diaiws have been made. Once each program has been compiled, 
Che Cask must be budlt before being executed. To do so, Che following cations 
must be in effect: ,'ASN CKx:*CI7:' where x is Che disk that contains Che 

CLOUD routines. Then, the coauoand: 'TKB @CL0UD' must be issued. This 

begins the task-buildi^ process. If a message 'Nor contigous disk space' is 
received, programs must be either eliminated or Purged &om the program disk 
In order to make exxsugh room for the task. 

3. Programming changes in the (XCUD routines: 

1. An addition of ve rsio n 7 constants (CONST?. FIN) . 

The subroutine CONST?. FIN is an exact duplicate of Che sxjbrouCine 
0CNST6.FIN except for the data statements. 

2. System library: 

The file LIB.CLB is us ed a s the system library during the task-building 
procedure. Once Che CONST?. FIN and KETIN.FIN routines were complete, 
these routines had to be addecT and replaced to the system library. 

a. Addition of CONST? Co the system library. 

To add CONST?, the command: 

'LBR LIB.OLB-<rNST?/lN' 
was issued. 

b. Replacement of KEYIN in the system library: 

To r^lace KEYIN, Che command: 

'LBR LIB.OLB^^CEYIN/RP' 
was issxied. 

3. CALL rcutlnes to properly use and assign the version 6 and ? constants. 
In order Co get the constants to work properly, two statements had to be 
added to the program PREPOS.FIN of the system: 

IF (NVERSN.BQ.6) CALL 00NST6 
IF (NVBRSN.EQ.?) CALL CONST? 

these two statements were added to PEIEPOS.FIN at the beginning of 
the program inmediately after Che statement OPFILE' . 
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APFBNDIX B 
GBAHT PBOPOSAL 

i 


m-aouD Ha phediciions for 

SPACE SHUTTLE LAUNCHES 


ABSTRACT 

Th« priaazy objtcclve of cht woric proposed hsrt will be co develop the 
cspebilitlcs of Che NASA/HSFC Multilayer Diffusion Model (MDM) to obtain 
in^cloud predictions of HCl concentrations in the Space Shuttle ground launch 
cloud. This will include documenting the procedures for running the MIRi on 
a POP 1145 and establishing the effect on in-cloud HCl concentrations using 
parameters characteristic of the standard aeteorologies encountered at 
Cape Canaveral. This information will then be used to establish an appro- 
priate aircraft sampling pattern prior co the March 1981 Space Shuttle launch 
CO both obtain representative oeasuremencs of in-cloud HCl concentrations 
and CO aid in verification of siodel predictions. 


1 


/ 
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schedule AMD DESCRIFTIOM OF PEOPOSED WORK 

IL« Approach uAod co carry ouc eha abova objacc;lvt would lacluda cha 
following phaaaa: 

PHaSE I: (Oaa ooach full tlaa ac JSC starting Aug. 1, 1980) 

a. During this tlaa it will ba nacassary to bacoma raacquaintsd with 
tha oparation of cha MDH as ir. is prograanad to run on cha PDP 1145 
and with cha modifications which hava baan mada by Larry Ray. 

b . Soma tiaa will ba davoctd co salaccing cha macsorological daca 
which will ba usad for in-cloud ccncantracioa pradiccions. Tha 
daca sacs should includa riprasantaciva macsorologias usad in 
previous studies (raf. 1) and from oora racanc and axcansiva 
casa scudias by Richard Silar. Daca should also Includa cha 
cast casa usad co varify cha modal in rafaranca 2. 

c. This daca will ba formacad (probably on disk) so chat niniaal 
afforc will bu raquiraa by JSC parsonnal whnn saccing up tha 
compucar co run frOTi a ramoca carminal. 

d. Tha faasibilicy of linking cha JSC compucar facilir/ co an 
incalligant carminal ac Kaarcey Scaca Collaga will hava baan 
cascad prior co assignmanc of this contract. Howavar, an important 
additional activity during this period will ba co further iniciace 
and check ouc all phases of che operation and input 7'd output 

for the on che PDP 1145 as accivaced by means of remote 
terminal from .Kearney State College. 
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PHASZ II. (Six Machs h clas tc Koartiay Scsct Collagt fcon S«pc. 1 co 
Fob. 28, 1980-81) 

4. Tha inlcl4l icsp co obC4inlag ehs In-cloud prtdlceloos will b« co 
dactcain* whachtc cha pocclon of ch« MlXl rMponsihl* for ch« 
pradiccioos is iacscc on chs program disk bslng ussd. 

b. If chs program La Incacc, Ic will b« oacassary co dscaralns chs 
approprlacs prograaml.ig co accsss cha oucpuc for various In-cloud 
lavtls. If cha subrouclnas for In-cloud concsn era cions of HCl 
ara ooc incacc, Ic will ba aacassary co dacarmina cha approprlaca 
progrsaming and ralncroduca ic. 

c. If case casa sacaorologlcal daca axlsc for which In-cloud 'rtdlccions 
hava pravlously baan oada, chay will ba usad in chis phasa co chack 
on cha raasonablanass of cha pradiccions onca chay ara obcainad. 
Anochar chack on cha accuracy of cha rasulcs taay ba co chack on cha 
consarvaclon of HCl ac various class afear cloud s cab llira cion. 

d. As soco as cha In-cloud pradiccions obcainad ara Judgad co ba 
raasonaola, savaral aacaorological casas characterizing diifaranc 
scabilicy ragiaans ac Capa Canaveral will ba used. The rasulcs 
froa chase casas will ba graphed co display cha vertical profiles 
of HCl as a function of claa and/or distance froa cha polnc of 
cloud scab lllza cions. Particular accancion will ba paid co the 
level ac which masi miia concantra cions of HCl occur for the differanc 
scabilicy classes. 

PHASE III. (Ihraa aumchs h claa ac Kaamay Scace College, March 1 - May 30, 
1981) 

Soma of the objactivas in Phasa II nay run over Inco Phasa III, because 
of the uncertainty in cha amount of time which will be required co successfully 
obtain in-cloud HCl predictions from cha MDM. 


3 
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a. As the proposed Space Shuttle Launch date In March 1981 approaches, 
the MDM predictions will be continuously updated as taore refined 
meteorological data relative to the conditions at the launch site 
are received. 

b. If the Space Shuttle Is launched on schedule and suitable HCl 
concentrations have been obtained, these results should be matched 
against the MDM predictions to verify the model for the meteorological 
conrlnlons existing at launch time. 

c. Crlterl.i for In-cloud air sampling patterns will be established for 
future launches based on a knowledge of the HCl predictions from 
Phase II and of the Cloud Stabilizations heights and its dependence 
on meteorological parameters as determined in previous studies . 

d. A final report covering all accivitias and results obtained over 
the contract period will be written. 
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